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We theoretically analyze a system where two electrons are trapped separately in two quantum 
dots on a suspended carbon nanotube (CNT), subject to external ac electric driving. An indirect 
mechanically-induced coupling of two distant single electron spins is induced by the interaction 
between the spins and the mechanical motion of the CNT. We show that a two-qubit iSWAP gate 
and arbitrary single-qubit gates can be obtained from the intrinsic spin-orbit coupling. Combining 
the iSWAP gate and single-qubit gates, maximally entangled states of two spins can be generated 
in a single step by varying the frequency and the strength of the external electric driving field. The 
spin-phonon coupling can be turned off by electrostatically shifting the electron wave function on 
the nanotube. 


I. INTRODUCTION 

As mechanical resonators with potentially high Q fac¬ 
tors and large resonance frequencie^iM^, ultra-clean single 
wall carbon nanotubes (CNTs) are promising systems for 
studying the coupling of the electron degrees of freedom 
to the mechanical motion of the resonatoil^lf^. On the 
other hand the properties of CNT such as valley degen¬ 
eracy and the curvature induced spin-orbit interaction 
attract much attentioiP^HUl two valleys in the elec¬ 
tron energy spectrum distiMuish semiconducting CNT 
from III -V semiconductors^. Qubits can be defined as 
the electron (hole) spins or the valleys in quantum dots 
(QDs) in CNTiSM] spin-orbit interaction due to the 
curvature of CNTs has been st udied both in theorjff^Ml 
and observed in the laboratorj^HIlI] The spin-orbit in¬ 
teraction pWs an important role as a source of spin 
decoherence^S and at the same time it allows the electri¬ 
cal control of the spin in bent CNT in a magnetic field as 
well as cooli ng of the CNT resonator using spin-polarized 
currentP2124| Furthermore, the coupling of the spin in a 
single QD and the deflection of the CNT was studiecP^, 
and the spin-phonon coupling, which is induced from the 
spin-orbit coupling where the tangent vector instanta¬ 
neously depends on the phonon displacement, provides 
a new platform for operating spins and quantized flex¬ 
ural mode^^. The read-out of the resonator vibration 
frequency and the detection of the single electron spin in 
the QD have been proposed based on the spin-phonon 

couplinJ^MH! 

Previously, we have proposed arbitrary single qubit 
gates using an electron spin in a single QD which lies in a 
suspended CNT making use of the spin-phonon coupling 
of the mechanical motion of the CNT^. For quantum in¬ 
formation and quantum computation, one-qubit and two- 
qubit gates are universaP^. Individual two-qubit gates 
that can form a universal set in combination with single¬ 
qubit gates are e.g. CNOT, VSWAP and iSWAI®®^. 
There are approaches related to inhomogeneous mag¬ 
netic fields to produce universal gates of spins and to 
achieve coupling of long distance spins in NV centers us¬ 
ing mechanical resonator^2lH2ll Compared with the use 



FIG. 1: (a) Schematic of the nanomechanical system where 
two quantum dots (QDs) lie in a doubly clamped, suspended 
CNT which is fixed by two supports at two ends. QDs are 
formed by the electronic potentials applied by the gate elec¬ 
trodes no. 1, 5 and 9 to form two QDs. Here, we assume that 
the fourth harmonic flexural mode of the CNT is excited. The 
standing wave in each QD is asymmetric and has one single 
electron trapped inside. The charged CNT is driven to vi¬ 
brate along the x axis by an external ac electric field applied 
by an antenna or the gate electrodes. A magnetic held is ap¬ 
plied along the z axis, (b) To switch off the coupling between 
two QDs, the QDs can be electrostatically shifted. The left 
QD between gates no. 1 and no. 3 as well as the right QD 
between gates no. 5 and no. 7 are both left-right symmetric. 


of inhomogeneous magnetic fields, electric fields are easier 
to control temporally and spatially. Universal quantum 
computation requires that arbitrary pairs of two qubits 
can interact with each other. It is usually not easy to ful¬ 
fill this requirement because long distance coupling can 
be very demanding. In the present paper we theoreti¬ 
cally study a two-qubit iSWAP gate and arbitrary single 
qubit gates in a nano-mechanical scheme where two elec¬ 
trons are trapped separately in two QDs on a suspended 
CNT. The indirect coupling of two distant single-electron 
spins in two separated dots is mediated by the vibra¬ 
tional motion of the CNT. A single-step preparation of 
maximally entangled states is obtained by combining the 
iSWAP gate and single-qubit gates. All quantum gates 
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proposed here can be controlled electrically. We show 
that the spin-phonon coupling in each QD can be turned 
off by electrostatically shifting the electron wave func¬ 
tion on the nanotube. We predict that arbitrary pairs of 
QDs can be coupled through the spin-phonon coupling 
on CNT with multiple QDs. When there are more than 
two QDs in the CNT, it is possible to couple arbitrary 
pairs of distant electron spins and at the same time turn 
off the spin-phonon coupling in the other QDs. 

This paper is organized as follows. In Sec. [IT] the 
nanomechanical system and the model Hamiltonian are 
introduced. In Sec. Illll a Schrieffer-Wolff transformation 
is applied to obtain the effective Hamiltonian in spin 
space and to obtain the iSWAP gate from the effective 
Hamiltonian. In Sec. |IV[ we determine the wave function 
of the qubit state by solving the Schrodinger equation 
and simulate numerically the fidelity of the maximally 
entangled states in a open quantum system by using a 
quantum master equation. In the Sec. we describe 
how to shift the electron wave function to turn off the 
spin-phonon coupling. 


II. MODEL 

We assume that a doubly clamped, suspended CNT is 
fixed on two supports at both end^^, see Fig. (a). 
Two QDs can be formed by applying proper voltages on 
gate electrodes below the suspended CNT. We further 
assume that a voltage is applied to the gates no. 1, 5 
and 9 so that the electrons are prevented from tunneling 
out of the CNT or from tunneling between QDs. The 
other gates can be used to tune the resonance frequency 
of the CNT. We assume that two electrons are trapped in 
two QD^, separately, and that an external longitudinal 
magnetic field By is applied along the z-axis of the CNT. 
An external ac electric field is applied by an antenna on 
the top or on the gates to excite the vibration of the 
charged CNT. 

Here, two single electron spins in two QDs are assumed 
to couple to the the vibrational motion simultaneously, 
hence these two spins are indirectly coupled via phonon 
exchange. We describe this system using the Hamiltonian 

H = Ho + Hu (1) 

= X! (2) 

i 

Hi = 2h\{a + a)) cos(uj t)+ a^) {a+i + a-i), 

i 

( 3 ) 

where z = 1,2 refers to two electrons in two separate 
QDs. The two spin states cross at the magnetic field 
B* « Aso/{2g,B) in the K valley of the ground state 
of a single electron QD in CNT^^SHUl^ where Ago is the 
spin-orbit coupling strength and /tb is the spin magnetic 
moment. We choose these two spin states as the qubit 


and assume that we are near the crossing point. The 
Zeeman splitting energy between qubits induced by the 
magnetic field is Hujqi = gg,B{Bi — B*). The quan¬ 
tized mechanical motion is described by the phonon mode 
with frequency uip and a (a^) is the phonon annihilation 
(creation) operator. We assume the system to be at low 
temperature hujp < k^T. Here is the Pauli z matrix 
of the electron spins and a±i are the corresponding spin 
raising and lowering operators. For simplicity, we only 
consider the fourth excited quantum harmonic flexural 
mode along the x axis of the CNT in the present paper. 

Due to the curvature caused by the vibrational mo¬ 
tion, the local tangent vector t of the CNT depends on 
the displacement coordinat^^, and it induces an inter¬ 
action between the mechanical motion and the electron 
spin. Spin-phonon coupling originates from the dynam¬ 
ical spin-orbit interaction cr ■ t(z) = cTz -|- {du/dz)crj;, 
where u(z) is the displacement at the coordinate point 
z, u{z) = + a'^), where /(z) is the waveform 

of the phonon mode and Iq is the zero-point displace¬ 
ment. For different QDs in the nanotube, the spin- 
phonon coupling strengths are gi = Ago {f {z))^lo/2\^. 
Here f {z) is the derivative of the waveform of the phonon 
mode. Considering the electron distribution, we obtain 
the average of the derivative of the waveform {f'{z))^ = 

(■^) where Dj(z) are the charge densities 
of two QDs and each quantum dot is between —li/2 
and li/2. We use realistic parameter L = 800 nm. 
Ago = 370 /reV, Iq = 2.5 pm and obtain the value of 
the coupling strength g/{2Tr) = 0.56 MHJ^. 

We assume that Di{z) are symmetric functions in the 
QDs. In this case, the spin-phonon coupling strength 
is non-zero when the parity of /(z) is odd in the QD. In 
other words, if the parity of the charge density function is 
even in a QD, to avoid canceling out the spin-orbit inter¬ 
action, there should be a left-right asymmetric standing 
wave in the qePsI. 


III. EFFECTIVE HAMILTONIAN FROM 
SCHRIEFFER-WOLFF TRANSFORMATION 

For a better understanding of the evolution of the 
spins, we derive the effective Hamiltonian in the sub¬ 
space of the spins. We assume that the difference between 
the phonon energy and the qubit energy is much larger 
than the spin-phonon coupling strength and the driv¬ 
ing strength, i.e., that ujp — ujqi ^ A. The Schrieffer- 
Wolff transformation can be applied when the subspaces 
with different phonon numbers are energetically well sep¬ 
arated. First, we obtain the effective Hamiltonian in 
the lowest subspace with zero phonon, then we can use 
the same method to obtain the effective Hamiltonian in 
all subspaces. The effective Hamiltonian from the time- 
dependent Schrieffer-Wolff transformation in the lowest 
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FIG. 2: The energy-level diagram of the combined two-qubit 
and phonon system. \(ja'n) denotes a state, where cr (cr') rep¬ 
resents the first (second) spin state and n is the number of 
phonons. The dashed lines denote the coupling strength of 
the external ac electric field A and the spin-phonon coupling 
strength g. The ac electric field is detuned from the phonon 
frequency by A. The coherent coupling between states I),), 0) 
and Iti 0 ) is mediated by the state I),), 1 ) through the driv¬ 
ing and the spin-phonon coupling. The effective single spin 
resonance coupling strength (5 and the effective two-spin in¬ 
teraction strength 7 are obtained by deriving the effective 
Hamiltonian for the n = 0 subspace and thereby eliminat¬ 
ing the \(ya'l) states with a Schrieffer-Wolff transformation 
(similarly for n > 0). Here F is the damping rate of the CNT. 


phonon subspace can be written a^^SEZl 


where 


HeS — 


II 

j Ci^zi x\^x2 

Hls = 

2 cos UjtJ2i 

fllOqi 

{2n -b l)hwqig1 


hXgiU}p{uj'^ - 2 ujI + 
%iff2Wp(-2a;2 -b ujI^ + ujI^) 

(u;2-a;2,)(a.2-a;2^) 


(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 


It is worth pointing out that in Eq. @ there is not only 
the coupling term which denotes coupling of two spins, 
but also the single electron spin rotation terms ax and 

To get rid of the time dependence in H^g, we trans¬ 
form the Eq.(|^ into the rotating frame with frequency 
w, using the transformation H^g = UHesW —iUU^ with 
JJ — gi(w/ 2 )t ^ assume Wp ~ Wq ~ W, A = Wp — w 
and A ^ g. The fast oscillating terms with can 

be dropped in the rotating-wave approximation. We ex¬ 
tend our analysis to the full phonon space and obtain the 
effective Hamiltonian from the Schrieffer-Wolff transfor¬ 
mation in the rotating frame (see Appendix [A|) 


O^iO'zi + Piaxi) + 7((T+i(T_2 + (T_iCT+2)) (10) 

i 



t[|JS] 


FIG. 3: The time evolution of the concurrence (color scale) 
as a function of the driving strength A with different initial 
states (a) |t4^) and (b) |tt)- The dashed lines denote the time 
evolutions of the concurrence with a fixed driving strength 
which are shown in Fig. When A = 0, in (a) the period 
between high concurrence peaks is to = ^(Wp — upp/(Ag^u-p), 
which depends on the strength g of the spin phonon coupling, 
while in (b), the spin-phonon coupling does not have any ef¬ 
fect on initial state |tt)- When X ^ g, the period of the 
concurrence peaks depends on the initial state and the peri¬ 
ods are t ~ (4/3)to in (a) and t « 4to in (b). The system 
parameters are chosen to be i2;p/(27r) = 1500 MHz, ij/{2'K) = 
tjq/(27r) = 1450 MHz, g/{2'K) = 0.56 MHz and n = 0. 


where 


CXi — 



(2n -b l)nujg,g'^ 




— UJ 


2 

qi 


huj 

T’ 


(11) 


and n = 0^0 is the phonon number operator. The energy- 
level spectrum is shown in Fig. Arbitrary single-qubit 
gates of the single electron spin can be obtained by com- 
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billing rotations about the a;-axis and the z-axi^^. The 
rotations about the z-axis of each QD can be adjusted by 
changing the driving frequency, and can be switched off 

by setting u) = a;qi(l — rotations about 

the cc-axis can be adjusted by choosing different strengths 
of the driving field, and it can be switched off by setting 
A = 0. From rotations about a:-axis (z-axis), X (Z) gate 
which is X = ax {Z = a^) can be obtained!^. 

We show how to obtain the iSWAP gate from the effec¬ 
tive Hamiltonian. We assume Wi = Wq, w = Wq and gi = 
g. The iSWAP gate is obtained in the absence of driv¬ 
ing, A = 0 where the phonon vacuum fluctuations cou¬ 
ple the two QDs. Choosing the appropriate pulse length 
t = n{ojp — ujq)/(4g'^ujp), we obtain the evolution opera¬ 
tor U' = in the basis {|tt), lU), lit), lit)} in 

the following forirPH 


A 


U' = 


V 0 


0 0 
0 i 
i 0 

0 0 e* 


\ 


/ 


( 12 ) 


We can see from Eq.(12) that the evolution operator is 


an iSWAP gate with relative phases. We can apply single 
qubit gates azi on the two QDs for t = 7 r( 2 /a;q — l/2a;p) 
to eliminate the relative phases between states |tt) and 

lit). 


IV. WAVE FUNCTION AND MAXIMALLY 
ENTANGLED STATES 

One can combine the iSWAP and single qubit gates to 
obtain any maximally entangled states from any initial 
product state. In our case, we can achieve a single-step 
preparation of maximally entangled state by adjusting 
the driving strength and driving frequency. The effective 
Hamiltonian in Eq.(10) is time independent. We solve the 

Schrodinger equation = H'^g{\n) Z> |V'(0)) 

in the subspace with phonon number n = 0 , where 
I^W) = Cl(t) Itt) + C 2 (t) lU) + C3(t) lit) + C4(t) lU) is 
the wave function of the qubit states with a initial prod¬ 
uct state. The exact maximally entangled states and the 
corresponding time points can both be obtained by solv¬ 
ing the Schrodinger equation. The Hamiltonian in the 
product basis is 


Kb = 


/ 2 a 13 P 0 

13 Q -i (3 

^ 7 0 /? 

\ 0 /3 /3 -2a. 


(13) 


We can use the four eigenvalues /ii(A, A,^) and four 
eigenstates |' 0 i) to obtain the general solution of the wave 
function |V;(t))) = Ai |V^i(t))) 

A 3 I'ljjsit))) -I-A 4 \'ip 4 (t))) where Ai is depen¬ 

dent on the initial state. 
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FIG. 4: The time evolution of the coefficients Ci{t) = 
{i\tp{t)} of the wave function and the concurrence with driving 
strength A/(27r) = 4 MHz for different initial states, (a) |t4): 
and (b) Itl")- (a^) The maximally entangled state obtained 
at ti = 26.5 gs is ~ (—0.09 + 0.23i) |tt) + (—0.13 — 
0.65i) |t4) + (-0.62 -t 0.23i) |;t) + (-0.16 + 0.21i) \U)- (b) 
The maximally entangled state obtained at t 2 = 76.8 gs is 
V>(t 2 ) « (0.52 -b 0.5H) in) + (-0.02i) |n) + (-0.02i) in) + 
(—0.48-b0.49i) in)- The other parameters are the same as in 
FiglU 


To quantify the generated entanglement, the concur¬ 
rence, as a measure of entanglement, is evaluated. The 
concurrence for pure two-qubit state can be written 
as C'(|V’(t))) = \ {'ilj{t)\ay iZ ay\ilj{ty) \ = 2 \ci{t)c 4 {t)- 
C 2 (t)c 3 (t)|. We show the concurrence as a function of 
time and driving strength for different initial states in 
Fig-i Since the concurrence reaches the value C = 1, we 
know that maximally entangled states can be obtained. 
From the solution of the Schrodinger equation, we obtain 
the time dependent C{t) of the concurrence. The pe¬ 
riod of the concurrence depends on the initial state, the 
driving strength and the spin-phonon coupling strength. 
Varying in time the driving strength, we can shift the 
maximum of the concurrence. In the following we take 
two initial states as examples in Fig. When A = 0, 
the period of the concurrence is t = 7 r(w^ — u:'^/{4g'^ujp) 
and only the iSWAP gate acts on the initial state. The 
iSWAP gate creates entanglement on the initial state 
Iti) in Fig. (a) but no effect with initial state |tt) 
^ (bjT When the driving strength is nonzero. 


in Fig. [| 













































5 


the electron spin resonances corresponding to X- gates 
are on. When 0 < A < 5 , the frequency of the electron 
spin resonances is slower than that of the iSWAP gate, 
therefore the period between high concurrence peaks de¬ 
pends on the spin resonance frequency. We can see from 
Fig. |3] that the frequencies of the blurred oblique lines 
in Fig. 1^ (a) and the bright oblique lines in Fig. i(b) 
are determined by the strength /3 of the single electron 
rotation Ux- When X > g, the frequency of the spin 
resonances is higher than for the iSWAP gate. The elec¬ 
tron spin resonances contribute a fast oscillation of the 
wave function. The periods of the concurrence depend 
on the strength of the iSWAP gate. Corresponding to 
the vertical strips in Fig. the periods of the concur¬ 
rence cycles are t « Tr{ujp — lo^)/(S g'^ujp) in Fig. [S (a) and 
t « -^uj^/ig^LOp) in Fig. ,0 (b). We also show the 
time evolution of the wave function with driving strength 
A/(27r) = 4 MHz in Fig. which corresponds to the 
dashed line in Fig. [^with X > g. Although the coeffi¬ 
cients of the wave function are fast oscillating due to the 
strong rotation ax, the envelopes of the time evolution of 
coefficients still correspond to the iSWAP gate. 

To include the damping of the CNT due to the cou¬ 
pling to a thermal bath at temperature T, we use a mas¬ 
ter equation for the nonunitrary dynamical simulation. 
The density of states of the other phonon modes is small, 
therefore we can expect a small spontaneous qubit relax¬ 
ation rate 1 /Ti and neglect it in the following. We obtain 
the master equation for the density matrix p, 


p = - Uh,p\ + {ub + l)r ( apa^ - -{a^a,p} 


+ nsF 


(14) 


where ub = _ 1 ) ig the Bose-Einstein oc¬ 

cupation factor, and T ^ g is the damping rate of the 
CNT. The phonons follow the Bose-Einstein statistics in 
the thermal equilibrium in the initial state of the density 
matrix that p = ^ 1 ^) (^1 ® 1 ^) (^ 1 ’ 

where Z = is the partition function. 

The total spin state is given by the partial trace over the 
phonons ps = Trphp. The fidelity relative to the two- 
qubit entangled state is defined as F = y /(vhl ps I'h). 

We solve the master equation to evaluate the fidelity 
of the entangled qubit state at a finite temperature in 
the presence of damping of the CNT. We choose the 
qubit state as |4') « (0.52 -|- 0.51*) |tt) + (—0.02*) |ti) -b 
(—0.02*) lit)+ (-0.48-1-0.50*) lit) which can be obtained 
at time tideai = 76.8 p,s with initial state |tt) at zero 
temperature with the parameters in Fig. In Fig. 

we plot the time evolution of the fidelity with fixed 
F = ojp/Q « 3x 10^ s“^ where Q « 30,000 is reachable in 
experimenlpHI. Because of the temperature T = 30 mK, 
the fidelity of the case with F = 0 is less than 1. From 
Fig.[5](b), we can see the peak of the fidelity with F = 0 is 
shifted from tjdeai = 76.8 ps. Although the iSWAP gate 
and the spin resonance mediated with virtual phonons 




FIG. 5: (a) The time evolution of the fidelity F — yj (4'|ps|4') 
with driving strength X/{2-n) = 4.0 MHz with a maximum 
near tideai = 76.8 ps obtained by solving the master equa¬ 
tion Eq. ( |14[ ) which takes into account the thermal equilibrium 
phonon bath at T = 30 mK and the damping rate of the res¬ 
onator. We magnify the circle in (b). (b) The fidelity for the 
case without damping is smaller than 1 due to the finite tem¬ 
perature. The shift between the peaks and the dashed line 
decreases when the damping rates increase. We truncate the 
phonon Hilbert space for n > 6. The other parameters are 
the same as in Fig. 


are not affected by the phonon numbers, the coefficient 
a of the rotation a^ at finite temperature is related to the 
phonons. The shift of the peak of the fidelity is larger 
at high temperature than at low temperature because 
the phonons follow Bose-Einstein statistics in the ther¬ 
mal bath (see Appendix [b|. 

We assume that the initial state is in the thermal equi¬ 
librium, so damping does not change the phonon dis¬ 
tribution. While a large damping rate slows down the 
transition process of the states with multiple phonons 
whose spin-phonon coupling strength of the states with 
phonon number n is y/ng, the fidelity increases with an 
increasing damping rate. In Fig. we plot the fidelity 
at tideai = 76.8 ps for |'I') with the initial state |tt) as 
a function of the damping rate. While the damping rate 
increases, the fidelity surprisingly displays a minimum 
when the damping rate approaches F « 2/3 and which 
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FIG. 6: (a) The fidelity F of obtaining the maximally en¬ 

tangled state |\[t) as a function of the damping rate F with 
initial state Itl")! taking into account the phonon bath in the 
thermal equilibrium at T = 30 mK. The initial state is Itl")- 
The stochastic resonance occurs at F — 2/3 = 0.56 x 10“® s“^. 
The other parameters are the same as in Fig. The fidelity 
for r = 0 is limited by the finite temperatnre T > 0. (b) 
The fidelity F of obtaining the maximally entangled states 
with different parameters as a function of the damping rate 
r. We separately modify the parameters A, A and g by 
a factor of 2, as \/{2n) = 2 MHz, g/2Tx — 0.28 MHz, or 
tj/27r = ujql2'K = 1400 MHz, while keeping other parameters 
as in (a). Hence the value of /3 « Ag/A is the same for all 
three cases, and half of its value in (a). The maximally entan¬ 
gled states in these three cases are different and obtained at a 
different time. The minimal fidelities in all these cases occurs 
at F = 0.28 X 10“® s“^, which shows that the maximum of 
the stochastic resonance is at F « 2/3 (vertical dashed line). 
We truncate the phonon Hilbert space for n > 6. 


we interpret as a stochastic resonanc^^. The damping F 
of the phonon leads to transitions between configurations 
with the same qubit states and different phonon numbers. 
The effective coupling /? is generated by the ac electric 
driving field with a large detuning as in Fig. From the 
effective coupling strength Ui in Eq.( 101, one can see that 
the phonon number n and therefore the effective coupling 
strength at fluctuate around their thermal average value 
with correlation time 1/F. These fluctuations do not sig¬ 
nificantly affect the coherent Rabi oscillation with Rabi 


frequency /3, except when their correlation time of the 
fluctuations is half of the period of the driving field. The 
stochastic resonance reaches its maximum at an optimal 
moderate value of the damping rate where F « 2/^. On 
the other hand the ac electric driving field is important 
for the electron spin resonance, but it does not have the 
effect of increasing phonon numbers for the large detun¬ 
ing A = ojp — oj. With the stochastic resonance of the 
driving field and the damping, the transition between 
states with the same spin and different phonon numbers 
is enhanced, which is detrimental for obtaining the elec¬ 
tron spin resonance and the ideal maximally entangled 
states. Therefore, we find a minimal fidelity at the max¬ 
imal stochastic resonance. 


V. COUPLING OF ARBITRARY QD PAIR IN A 
QD ARRAY 

Universal quantum computation requires that arbi¬ 
trary pairs of qubits can be coupled. We extend the 
case of two QDs to several QDs on the CNT, with one 
single electron trapped in each QD. To couple an arbi¬ 
trary pair of QDs, the coupling between different pairs of 
QDs should be controllable. In other words, we should 
be able to turn on and off the two qubit coupling be¬ 
tween any arbitrary two qubits on the CNT. The coupling 
of two qubits is bridged by the spin-phonon interaction, 
hence it is possible to cut the coupling by breaking the 
spin-phonon coupling. We show in the following that 
the interaction between any arbitrary two qubits can be 
switched off and on by controlling the spin-phonon cou¬ 
pling in each QD. 

The spin-spin coupling in two QDs is induced by the in¬ 
herent spin-phonon coupling in each QD. As we have dis¬ 
cussed in Sec. |llj under the precondition of a symmetric 
charge density function, the spin-phonon coupling is can¬ 
celed if there is a symmetric distribution of the phonon 
waveform in the QD. By adiabatically changing the volt¬ 
ages which form the QDs, we can tune the location of 
the QDs to lie at the anti-nodes of the vibrational stand¬ 
ing wave. When the phonon waveform is symmetric in 
the QD, the spin-phonon coupling is eliminated. In other 
words, we can turn off the spin-phonon coupling by elec¬ 
trostatically shifting the electron wave function on the 
CNT. The left dot, for example, is between gates no. 1 
and no. 3 in Fig. (b), therefore both distributions of 
the electron and the phonon waveform are symmetric in 
the QDs and the coupling strength of the spin-phonon 
coupling is zero for each QD. To switch the interaction 
on, we can electrostatically shift the electron wave func¬ 
tion to have an asymmetric phonon waveform in the QDs 
when charge density function is symmetric. Therefore, 
arbitrary pairs of QDs could be coupled from multiple 
QDs in CNT. 

It is possible to produce the X, Z, and iSWAP gates, 
separately or together, by adjusting the strength and fre¬ 
quency of the driving field and the positions of the two 
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QDs in CNT. A series of one-qubit gates and iSWAP where U{t) is a unitary transformation. If one writes 
gates is sufficient for arbitrary quantum computation. U{t) = where S{t) = oc 0{Hi), the trans¬ 

formed Hamiltonian at second order is 


VI. CONCLUSIONS 

In summary, we have studied a nanomechanical sys¬ 
tem consisting of a suspended CNT where two separated 
single-electron spins in two QDs are coupled indirectly 
via the exchange of virtual phonons. The CNT is driven 
by an ac electric field in a parallel static magnetic field. 
The indirect coupling of the two spins is provided by the 
simultaneous coupling between the two spins and the vi¬ 
brational mode of the CNT. We show that an iSWAP 
gate can be obtained by analyzing the effective Hamil¬ 
tonian derived from the time dependent Schrieffer-Wolff 
transformation and the time evolution operator when the 
driving electric field is off. Arbitrary single-qubit gates 
can be obtained in each QD by adjusting the strength and 
the frequency of the electric driving field. The iSWAP 
gate can be switched off when suppressing the spin- 
phonon coupling by electrostatically shifting the electron 
wave function on the CNT. Combining the iSWAP gate 
and single-qubit gates in the double QDs in the CNT, a 
universal set of quantum gates can be built and maxi¬ 
mally entangled states of two spins can be generated in a 
single step by varying the frequency and the strength of 
the external electric driving field. In this way, arbitrary 
pairs of distant spins in a QD array could be coupled. 
The fidelity for obtaining a maximally entangled state at 
a fixed time at finite temperature can be highly increased 
by increasing the damping rate of the CNT resonator. 
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Appendix A: Effective Hamiltonian from 
Schrieffer-Wolff transformation 

In Appendix we show how to obtain the effec¬ 
tive Hamiltonian from the Schrieffer-Wolff transforma¬ 
tion and obtain the form in the interaction picture. We 
can obtain the time-dependent Schrieffer-Wolf transfor¬ 
mation from the Schrodinger equation as 

H,s = UHU^ - iU{dtU^), (AI) 


i7eff = i7o + 0{H^) + 0{Hl). (A2) 

The first order terms 0{Hi) are eliminated to obtain a 
block-diagonal Hamiltonian that 


Hi + [S{t) ,Ho]+iSit) =0, (A3) 


from which the expression of S{t) can be obtained. Sub¬ 
stitute the expression of S{t) in Eq. (Al), we can obtain 
the transformed effective Hamiltonian Hes- For simplic¬ 
ity, we assume = gs = g, Wqi = Wq, A = Wp - w 
and A ^ g. We obtain the effective Hamiltonian in the 
lowest subspace with 0 phonon from the time-dependent 
Schrieffer-Wolff transformation in matrix form 




FIG. 7: Time evolution of the fidelity near tideai = 76.8 gs 
with driving strength X/{2ti) — 4.0 MHz at finite temperature 
T = 10 mK in (a) and T = 100 mK in (b). Inset: Fidelity 
over a larger time range; the circle denotes the area of the 
main plot. The other parameters are the same as in Fig. 

















































HeS = 


—^p+^q) 

2 A^cjp (uj^ ~2a;p +qjq) cos ujt 

(i^ 2 -qjg)(w 2 -u; 2 ) 

2A(;a;p (u;^ —2tJp +‘^q) cos ut 
2 g^qjp 


2 A 5 a;p(a;^ — 2ct;p+ci;q) cos ujt 2Xgujp{uj^ — cos tut 


(q; 2 -q; 2 )(qj 2 _qj 2 ) 

0 

2 g^qjp 

‘^p-‘^q 

2 A 5 a;p(a;^—2u;p+u;q) cosujt 

(q; 2 -p; 2 )(q; 2 _q, 2 ) 


(q;2-q;2)(q;2_q;2) 

2 g^qjp 

“p-‘-Jq 

0 

2Agu;p(a;^ —2£A;p+£A;q) cos 


2 g^ccp 


2Apa;p (cj^ —2cL)p+CL)q) cos cji 

(w 2 -qj 2 )(qj 2 _qj 2 ) 

2 A( 7 £A;p(a;^ —2a;p+a;q) cosujt 
(w 2 -w|)(qj|-qj 2 ) 
q(2g —tCp+qJq) 


and 


CUp ^^yu/pV^LU i,u/p-ru/q^ u/i. i.y\yujp\^LU uj l u/pTLUqy 

\ (q;2-q;2)(q;2_q,2) (q;2-q;2) (^2 _q,2) -CC^+CC^ / 

We transform Eq. (A4) into the interaction picture with respect to Hq. The fast oscillating terms with e=*=h‘^+“q)t 
± 2 iui^t i^g dropped in the rotating-wave approximation. We find 


His/^ =( + ^Bz) 

-w 2 + a ;2 

Xgujpiuj'^ - 2 cj 2 +a; 2 ) 
(a ;2 - w 2 j(a ;2 - w 2 ) 

2 g'^U] 


(e“* + e-“*)(CTA+e“‘>* + aB+e“<>* + aA-C-^'’* + 


(A5) 




(cta+ctb- + 0-A-crB+ + aA+CTB+e + (JA-aB-e 2*‘^q*), 


We obtain the effective Hamiltonian in the rotating frame 
with U = e*“*: 

His/n = - ~ i^Az + aBz) 


(2wp - A)A 
Xg2ujp 
(2wp - A)A 

2 g^cup 

(2wp - A)A 


i^Ax ~t” ^Bx) 


(cta+CB- + CrA-CTB+). 


(A6) 


the fidelity without damping and the ideal time point 
Tideai = 76.8 fjLS increase while the temperature increases. 
We can compare the two cases at T = 10 mK and 
T = 100 mK in Fig. The shifting is obvious at 
T = 100 mK but very small at T = 10 mK. This is be¬ 
cause the phonons follow the Bose-Einstein distribution 
in the thermal bath and the coefficient a,- of the rotation 


az in the effective Hamiltonian in Eq. (10) is related to 
the phonons. 


Appendix B: Fidelity with T = 10 mK and T = 100 mK 

When we consider the case with the phonon bath and 
the damping effect, the detunings between the peak of 
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